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Abstract 

We investigate the Exact Renormalization Group (ERG) description of (Z2 invari- 
ant) one-component scalar field theory, in the approximation in which all momentum 
dependence is discarded in the effective vertices. In this context we show how one 
can perform a systematic search for non-perturbative continuum limits without 
making any assumption about the form of the lagrangian. The approximation is 
seen to be a good one, both qualitatively and quantitatively. We then consider 
the further approximation of truncating the lagrangian to polynomial in the field 
dependence. Concentrating on the non-perturbative three dimensional Wilson fixed 
point, we show that the sequence of truncations n = 2, 3, — , obtained by expanding 
about the field (p = and discarding all powers (p 2n+2 and higher, yields solutions 
that at first converge to the answer obtained without truncation, but then cease 
to further converge beyond a certain point. Within the sequence of truncations, 
no completely reliable method exists to reject the many spurious solutions that are 
also generated. These properties are explained in terms of the analytic behaviour 
of the untruncated solutions - which we describe in some detail. 
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In ref.[3] we claimed that a sequence of truncations of the field dependence of 
the ERG[1][2] do not work in general (in the sense of providing dependable and 
in principle arbitrarily accurate results), because these do not converge beyond a 
certain point, and because no completely reliable method exists to reject the many 
spurious solutions that are also generated. In this letter we verify this claim for the 
relatively simple O(p ) case described in the abstract. We compute the truncations 
to high order (n = 25), and show how this behaviour may be understood, accurately 
and at a deeper level, in terms of the analytic behaviour of the untruncated solutions 
- which is described here in much more detail than was possible in ref.[3]. (We 
must emphasise here the distinction between a momentum/derivative expansion in 
which higher space-time derivative terms are discarded but no approximation is 
made in the field dependence - these approximations do appear to converge [4] [3] 
and truncations of the field dependence of the lagrangian which in general do not 
converge and can give even qualitatively wrong results) . 

We point out that the analytic properties of the untruncated solutions allow for the 
possibility of searching, within the momentum/derivative expansion [3] -[5], system- 
atically through the infinite dimensional space of non-perturbative lagrangians for 
new continuum limits. It need hardly be stated that, firstly very little is known 
about the possible existence of continuum theories in such a space, and secondly, if 
new fixed points were found, they could have profound implications. We have done 
such a search for O(N) scalar field theory in D = 4 and 3 dimensions, for the cases 
N = 1, 2, 3, 4, in the 0(p ) approximation. However, we found there to be only the 
known fixed points at this level: Gaussian for D = 3 and 4, and the Wilson fixed 
points in 3 dimensions. 

Finally we construct two better methods of approximation by expansion. The 
simplest way to solve the O(p ) equations, is directly numerically [3] however 
(see also later, and ref.[10]). The truncations, if carefully interpreted, can give 
moderately accurate results - thus the simpler low order truncations may be of some 
use in situations where more reliable and accurate calculations are prohibitively 
difficult to perform. This situation seems very reminiscent of approximations (also 
involving truncations of the operator basis) to the "real space renormalization 
group" investigated in the late 1970's[12]. For recent work on approximations to the 
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ERG see for example refs.[3]-[ll]; From Alford's work [8], and the numerical solution, 
it is clear that expansion of <p around the semi-classical minimum of the effective 
potential, results at low orders of truncation, in convergence to three decimal places 
or more for the O(p ) v (or the other exponents related by scaling). It would 
be interesting to better understand the reliability of this method and its limiting 
accuracy (also for uj). Its behaviour no doubt is otherwise similar to the truncations 
we discuss, and for the same underlying causes. Effectively this expansion was part 
of the calculation in ref. [9] , where the authors claim to obtain quite accurate values 
for exponents of 3D iV-vector models. At higher orders in the momentum expansion 
we expect that truncations become more limited in accuracy and reliability, if only 
because it ought now to involve the much more complicated calculation of the 
expansion of several functions[3] in polynomials, each with their own limitations 
in accuracy. Indeed, even in ref. [5], where we will consider truncations in the 
field dependence of higher momentum terms (0(p n ) with n > 0) only, making no 
expansion in the potential, the results support this hypothesis. 

Here we will concentrate on the case of sharp cutoff and O(p ), i.e. the remarkable 
equation[13] first studied without further approximation, by Hasenfratz and Hasen- 
fratz[10]. We start with a simple alternative derivation. 1 The equivalent smooth 
cutoff equation [3] is qualitatively very similar so that much the same analysis, and 
all our general conclusions, apply equally well to the smooth case. We work in D 
euclidean dimensions with a single real scalar field (p. The partition function is 
defined as 



Aa = A(q, A) = 6 e (q,A)/q 2 is a free massless propagator times a smooth (every- 
where positive) infrared cutoff function e {q, A). The sharp cutoff limit is given by 
the Heaviside function: 




(1) 



e e (q,A)^e( q -A) 
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From (1) we derive the flow equation for Wa- 




For more information on this, see ref [4]. 
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which on rewriting in terms of the interaction part of the Legendre effective action 
via T\[(p] + ^(p.A^.cp = -W A [J] + J.<f, <fi = 5W A /SJ, gives 



as™ 
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"2 tr 
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Ta generates the one particle irreducible parts of the Wilson effective action 
(Wegner-Houghton effective action[2] in the sharp cutoff limit) [4]. If we now make 
the approximation of discarding all momentum dependence from Ta in (3), we 
obtain 



d_ 
OA 



d$ £ (q, A) 
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OA e e (q,A)[l + e e (q,A)V»(<p,A)/q*] 



(4) 



where we have introduced the potential through T\ = Jd D x V(y(x), A). Primes 
denote differentiation with respect to (p. The integrand (in curly brackets) in (4) 
may be written 



-|- {in e £ (q, A) - In [l + 9 £ (q, A)V" A)/q 2 } } 



dA 



The first term above yields a field independent vacuum energy which we drop from 
V. Taking the sharp cutoff limit (2) we thus obtain for the integrand 

8(q-A)\n(l + V"( l p,A)/q 2 ) . 

The integral in (4) is now trivial. Factoring out the scale A by writing ip i— > 
A D / 2 "V/C ^(C"VA D/2_1 ,A) i-> (- 2 A D V(v,t), and t = ln(A /A), where the 



factor C = (4n) D / 4 ^jT(D/2) is chosen for convenience, we obtain the advertised 
equation[13] [10]: 

^V(<p, t) + (D/2 - \)<pV\tp, t) - DV(cp, t) = In [1 + V'\<p, t)} . (5) 

The anomalous dimension r\ = in this case since a non-zero rj results from non- 
trivial wavefunction renormalization. The reader may be puzzled as to why we 
obtain the same equation as that for the Wilson effective potential[10] given that 
ours is the Legendre effective potential. In fact these are one and the same, since at 
zero momentum all vertices of the Wilson effective action Sa coincide with Ta [4]. 
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Now we set D = 3. From (5), a fixed point effective potential V((p, t) = V(ip) must 
satisfy 

^y'M-3y(^) = in[i + y»] , 

and V"(0) = 

(by </? <-> — invariance). At first sight these equations appear to have many 
solutions, which may be parametrized by V(0) = — ln(l + cr)/3, a A 2 being the 
semi-classical effective mass-squared. Actually this is not the case [10], because all 
but two solutions end at a singularity of the form V(<p) ~ 2(1 — <p/<p c ) hi^c — ip), <p c 
some positive constant, or more precisely, as a series in decreasingly singular terms, 

-r j. "\ / \ ( 3 2 25 3 5y?2 3169 4 \ 

V = HX) [ X ~ 8 X - 432* - W4 X + 27648* + ' ' ' J (?) 

+ ln(x) 2 (^x 3 - ^ -■••)+ 0(ln(x) 4 x 5 ) , 

where x = 1 — ip 2 We now justify this statement, by dividing the behaviour of 
y(</?) into three classes. 

First of all, if V ends at a singularity, it, or some derivative of it diverges there. By 
considering how the various terms can balance in (6), one sees that a singularity 
must be of the form above. Secondly, if V does not end at a singularity but instead 
it and its first two derivatives tend to a limit (±oo included) as <p — > oo, then again, 
considering the balance of terms in (6) one sees that either V(ip) — > or V((f) 
satisfies 

V(<p) = A^-Un^- 2 -- 1 - ln(30A) - + 0(l/<p 6 ) , (8) 

for some positive constant A, as <p — > oo. Linearizing in (6) about the first possibility, 
i.e. setting V(<p) h- > V^(<^) +5V(y?), one finds that 5V must behave for large ip as 
a linear combination of </? 6 and exp(</? 2 /4). Since both of these corrections are 
excluded if we require also V((p) + 8V((p) — > as — > oo, we conclude that V(^) is 
identically zero in this case, which indeed is the trivial Gaussian solution to (6) . (For 
a study of perturbations about the Gaussian in eqn.(5), see ref.[10]). Linearizing 
about the second possibility, i.e. taking V((p) to be as in (8), one finds that 8V 
must behave for large ip as a linear combination of ip 6 and exp(5A</? 6 /2). Now the 
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first correction merely perturbs the coefficient A, while the second is excluded if 
V + SV also satisfies (8). It follows that the space of solutions obeying (8) divide 
into isolated one-parameter subsets, each parametrized by A. For both possibilities 
the exponentially growing corrections were the result of linearizing the singular 
behaviour (7), since they involve balancing the same terms in (6). Thirdly, V or 
one of its first two derivatives may be defined for all finite real <p but not tend 
to a limit as tp — > oo. Studying (6) one sees that this requires solutions to become 
infinitely oscillatory as ip — > oo. We do not supply a proof that this does not happen, 
but we feel confident we can rule it out because we saw no hint of it in our extensive 
numerical and analytic investigations (see below). 




0.0 1.0 2.0 3.0 





Fig.l. The solution V(p) which approximates the Wilson fixed point. 

Therefore, apart from the trivial solution V((p) = 0, any global fixed point solution 
of (6) must satisfy two boundary conditions: (8) and V'(0) = 0. We thus expect at 
most a countable number of such solutions; we find only one. It may be characterized 
by o = a* = -.46153372 • • • (or A = A* = .003033 • • •) and is displayed in fig.l. (We 
describe the search below). This is the same solution as in ref.[10] of course and is a 
fair approximation to the Wilson fixed point. To obtain the critical exponents and 
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operator spectrum one linearizes about this solution in (5), i.e. one writes V(cp, t) = 
V((p) + 5V((p, t), with 5V((p, t) oc v((p) exp(At). If one excludes the exponentially 
growing solution (which is exp(5A</? 6 /2) again) one obtains a discrete spectrum, 
and deduces from the eigenvalues the critical exponents, in reasonable agreement 
with the best estimates, as already described in ref.[3] (cf. also [10] and later). It 
is surely the case that the exponentially growing solution is again a linearization of 
singular behaviour, this time in the flowing solution V(ip, t). 

Note that it is wrong to conclude from fig.l, e.g. by naive semi-classical considera- 
tions, that the Wilson fixed point describes a spontaneously broken theory: indeed 
this is inconsistent with the fact that the field theory is scale invariant at this point. 
The remaining quantum corrections for momenta q < A, which generically give 
positive contributions to the mass-squared, in this case exactly cancel the negative 
mass-squared a* A 2 , and as A — > one recovers the complete Legendre effective 
potential as V((p) = A*(p 6 (from (8) in unsealed units). Note also that it is not 
necessary to consider separately the question of the physical stability of solutions 
to (6). This is because there are no solutions that are unbounded from below 
(contrary to the statement in ref.[10]), while solutions that end on the singularity (7) 
are clearly physically unacceptable: the potential does not diverge at <p c , it simply 
ceases to exist - or is complex thereafter. (In the Gaussian case the stability is seen 
to hold automatically once perturbations about the fixed point are considered, as a 
consequence of the In in (5)). 

The general structure of the solutions is as follows. For a close but less than a*, 
the solutions look very similar to fig.l but end at some cp = <p c in a singularity (7), 
as is most easily seen by plotting V" ~ 2/(p c /((p c — (f). As a approaches cr* from 
below, <p c moves out along the real axis, but very slowly, so that for <p c > 3 we 
require a* — a < A where A « .005, while a must approximate a* to high precision 
for say cp c <i 4. For a at the same proximity but above a* the singularity splits 
into a complex conjugate pair close to the real axis, with real positions Re(</? c ) 
in approximately the same place. These move closer to, and out along, the real 
axis as a — > a+. The distance from the real axis is also a sensitive function of 
a — a*. (Perhaps the position of these singularities is given by <p c ~ (a* — cr) _T , 
when a ~ a*, cr^a*, for some small positive constant exponent r). On the real 
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axis, at cp ss Re(<^ c ), V" turns over steeply and drops rapidly to a value just greater 
than —1 and approximately constant over a large range. In this range V slowly 
"rolls over" as 

1 1 ff 5 

V(ip) = -- InS + cip- -if 2 + 5 / dip (<p - ip) exp{V> 2 - -ciJj} H (9) 

3 2 Jo 2 

where <p> = c + 0(6) is the point inside the range where the potential reaches a 
maximum, 5 is small (5 — > as a — > cr+), and |<p — c| « \J\n(l/5). Eventually, at 
some position cp = <p' c outside this range, V encounters another singularity of the 
form (7). 

For a = —1 + 5, 5 — > + , the solution V(<p) is governed by (9) with c = 0, ending in 
a singularity cp' c ~ ^ln(l/<5). As a — > 0~, the singularities move out to infinity in 
such a way that V(<p) tends pointwise to zero. For a > 0, V(y>) grows monotonically 
(with real increasing cp) ending in a singularity (7), which moves out slowly as a is 
increased. 

Returning to the true solution, V((p) at a = a*, we note that V has a four- fold 
symmetry in the complex plane: complex conjugation x (ip <-> — <p). Factoring out 
this symmetry, if one carefully integrates out along rays <p = r e^, with < $ < 7r/2, 
one can determine the position of the closest singularity. It appears at r = r* = 3.12, 
d = t?* = .2577T. (There are others with r > r*). 

We see that it is possible to make a systematic search for new continuum limits 
without making any assumption about the form of the bare potential. Indeed if such 
a potential can be tuned to a critical point where a continuum limit is recovered, 
then the corresponding Wilson effective potential must satisfy (6) at that point. 
On the the other hand, if we find an acceptable solution of (6), then because such 
a solution is scale invariant, this equally well serves as the critical bare potential. 
Thus a search over the infinite dimensional space of bare potentials reduces to a 
one dimensional search for effective potentials obeying (6) that do not end in a 
singularity. In fig. 2 we plot the inverse of the position of the real singularity against 
o. The Gaussian and Wilson fixed points are clearly seen as sharp downward spikes, 
at o~ = and a = a* respectively. We have performed the equivalent search in O(N) 
scalar field theory in D = 3, 4 dimensions for iV from 1 to 4, as mentioned in the 
introduction - finding only the expected fixed points. (The relevant equation for 

7 



general N is given for example in ref.[10]). Of course the restriction to considering 
only general potentials is a result of our approximation; at higher orders in the 
momentum expansion a larger space of lagrangians can be searched. 

0.50 
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a 

Fig. 2. The inverse of the position of the real singularity (p c as a function of a, 
sampled with a stepsize 5a = 1/120. The downward spikes reach further towards 
zero with a finer mesh. The various features are explained in the description of the 
analytic structure given earlier. 

(Note that equation (6) is stiff [14], the higher order equations[3] more severely so. 
To obtain an accurate representation of the solution at a = a* one can binary chop 
between the slow rollover behaviour (9) for a > a* and the singular behaviour (7) 
for a < a*, but a more efficient method is to require (8) at some large value of 
ip. One can then either shoot to the origin 2 - determining A to zero V'(0), or use 
relaxation.) 

We can now give an intuitive explanation, based on the simplified context of (6), 
for why the truncations in fact converge at first but then cease to further converge 
beyond a certain maximum n. The point is that, as well as the true solutions, 

2 For N 1 one would need to shoot to an intermediate fitting point. 
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there are many 'bad' solutions with singular field dependence on the real axis. 
Very bad solutions have (real or complex) singular field dependence very close 
to the origin (p = 0, causing the coefficients of (p m (m-point Green functions or 
vertices in general) to diverge rapidly with m according to the appropriate radius 
of convergence. Naturally, the polynomial field dependence of the truncations, for 
which the 2n + 2 vertex vanishes, 3 tend therefore to better approximate the Taylor 
series of a true solution. Increasing n will tend at first to further improve the 
approximation, by in effect ensuring that the singularities are forced further from 
the origin. However a non-trivial true solution also has singularities for complex ip 
at (and in general beyond) some radius \ip\ = r*. Therefore the truncations cannot 
be expected to converge to better results than would be obtained from 'moderately 
bad' solutions with singular field dependence only at or beyond the radius r*. 

Now let us make this argument much more precise. If we wish to ensure that the 
potential V(cp) is an even function, then the Taylor expansion must be done about 
the origin. It is helpful to write 

n ,) = 4,„( 1 + ff ) + i V + 4E^^ • (10) 

Plugging this into eqn.(6) we deduce «4 = — + a), a 6 = ^a(l + a)(l + 7a), 
a 8 = —jgcr 2 (l + cr)(l + 3a), • • •, and for k > 4 that a 2 k is given as a 2 times (1 + a) 
times a polynomial in a - as follows from D = 3 being an upper critical dimension[7], 
the slow rollover behaviour as a — > — 1, and the recurrence relation 

«2fc+2 k - 3 A (-4) m ~ 1 ^ a 2 fc x +2 • • • «2fc m +2 

l + a-2k(2k-lf 2k + £- 2 m ^..^ (1 + a)- 

fci H h&m — & 

respectively. The n th truncation is defined by setting «2n+2(o") = 0. We concentrate 
on large n, where the solutions a that result from this, can be understood from the 
asymptotic expressions for a 2 k- To leading order in 1/k, if the closest singularities 
(7) are on the real axis at (p c = ±r one has ct 2 fc ~ l/r 2fc , while if they are complex 
then there are four in the form <p c = ±r exp(±i^) and one has a 2 k ~ 2 cos(2ki}) / r 2k . 
(In fact the a 2 k asymptote to these expressions even for quite small k; for example 



For a discussion of truncations in general see for example ref[4]. 
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one finds for a = a*, that the asymptotic expressions give the right sign for k > 2 
and are within a factor 2 for k > 7). Recalling the analytic behaviour of V(ip) as 
a function of a ~ a*, we see that for a < a* — A the coefficients ct2n+2(c r ) are all 
positive, so there are no solutions a in this region for large n. On the other hand 
for a = cr* the coefficients are controlled by the singularities <p c = ±r* exp(±i#*) 
and thus the signs of the 0^+2(0*) are very closely four-fold periodic in n in the 

pattern + H , as a consequence of ~ 7r/4. For the negative pair, it follows by 

continuity in a that 0:271+2(0") must vanish for some a in the range a* — A < a < a* . 
For the positive pair, one has to look at a > a* + A; here the closest singularities are 
complex with an angle $ which is a rapidly increasing function of a. This implies 
that there is always a zero of 02,7+2(0") close to but greater than cr* + A. Together, 
these solutions a are the real ones that best approximate the Wilson fixed point. 
They are shown in fig. 3 up to n = 25. (For the higher n one must work to an 
accuracy of at least 20 significant figures to avoid round off errors). One clearly 
sees the four-fold periodicity with amplitude ~ A, and that, as expected, the upper 
solutions are slightly worse approximations. Indeed the average of the last four 
solutions gives a = —.4607 which is greater than cr* by A/6. 

The four-fold periodicity is transferred to the critical exponents, which may be 
computed by linearizing, in (5), as 02fc(t) = a2fc(o") +£02k (where we have used 
separation of variables). The neatest method of computing the eigenvalues A is by 
the linear recurrence relations between the /?2fc) imposing (32 n +2 = 0. All truncations, 
but n = 23, have one positive eigenvalue as required, which yieldsfl] v = 1/A. The 
results are shown in fig. 4. One clearly sees again the limiting periodic behaviour 
about the exact solution [v = .6895), this time with amplitude .008. The exponent 
for the first correction to scaling is given by 00 = —A where A is the least negative 
eigenvalue. This depends much more sensitively on the approximations and is even 
complex for n = 19, 22, 23. It bounces about the exact answer (ui = .5952) with 
an amplitude ^ .15. These truncations have been considered before, in refs.[6][7] 
to n = 11, 8 respectively, but without the deeper understanding it was possible at 
these orders to interpret the numerical results as indicating convergence. (The much 
worse behaviour for O(N) scalar field theory with N = 2, 3 [7] is surely due to the 
exact solution having complex singularities much closer to the origin). 
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Fig. 3. The solutions a that best approximate the exact answer (a = a*, shown 
as a continuous horizontal line) for the truncations n up to n = 25. The solutions 
not displayed (n = 2, 3) lie outside the range on the a axis. The dotted lines are 
a = a* ± .005. Recall that A « .005. 

Of course these approximate solutions are not the only solutions for a. For n too 
small the asymptotic pattern has not set in. Even for n large, there are many 
'spurious' solutions «2n+2 = in the range a* < a < 0, which must be there by 
continuity arguments; away from the boundaries of this range, one finds that they 
slowly drift leftwards with increasing n (asymptoting towards a = cr*), reflecting 
the fact that the closest singularities have an angle which is a slowly increasing 
function of a in this range. Looking only at the truncations, how can one tell 
that these solutions are spurious? There is surely no completely reliable method. 
There is no good reason to require the truncated potential to be bounded below and 
indeed the first cases to violate that criterion are truncations n = 6, 7 in figs. 3, 4; 
nor is it compelling to assume the solutions must be real: the approximations are 
bad (compared to their neighbours) for truncations 22 and 23, unless one chooses 
certain solutions with small imaginary parts (a = —.4572 + .0059z and —.4566 + 
.0027i respectively - only the real parts are shown in figs. 3, 4). For n = 23 even 
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Fig.4. The exponent v for truncations up to n = 25. The exact answer is shown as 
a horizontal line. 



the requirement that the approximation should have only one relevant direction, 
breaks down: there are no such solutions. The case we chose has A = 1.472 + .0253z, 
which we use to compute v = .6791 — .0117z, and a less relevant direction with 
A = .509 + 1.695z. The spurious solutions generally, but not always, have about the 
same number of positive eigenvalues A as negative eigenvalues: which is many for 
large n. The best one can do [6] to eliminate spurious solutions, is to look numerically 
for convergence/stability with increasing n. That this is not good enough is nicely 
demonstrated in ref.[6] where the slow drift of a sequence of spurious solutions a, 
with two relevant eigenvalues, is mistaken for convergence to a tricritical point. 

Finally we mention two better expansion methods. The first is the analytic 
equivalent of shooting to an intermediate fitting point, namely we require the Taylor 
expansion (10) and its derivative to agree with the asymptotic expansion (8) and its 
derivative, at some given intermediate point cpf. In this way one obtains just one 
solution as expected, with a ~ a* and A m A*. By comparing a pair of truncated 
Taylor expansions with a pair of truncated asymptotic expansions over a range of 
(pf (to, at the same time, estimate the error and determine the 'best' (pf) one can 
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extract moderately accurate bounding values for a*, for example within a range 
of width .04 by using Taylor series to (p 8 and (p 10 and asymptotic series to (p~ 10 
and </?~ 14 . This method will be reliable providing the asymptotic series is accurate 
for ip > a (when A ?z A*) for some a such that a < r*, which is certainly case 
here. Presumably it works in general if there are no singularities in the region 
a < Re (<£>)< oo, except that it cannot give unlimited accuracy since the expansion 
(8) converges only in the Poincare sense. 
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Fig.5. Results for the exponent v against n, in an expansion method that utilises 
the fact that i?* 7r/4. The exact answer is shown as a horizontal line. Note the 
much finer vertical scale compared to fig. 4. 

In the second method we assume that a* has been determined to the accuracy 
required of the eigenvalues. Linearising about the fixed point position of the nearest 
singularity as ft = ft* + eipe^ and r = r* + ese^, one obtains the asymptotic 
behaviour of (32k- 

P 2k ~ -4k {s cos(2kft*)/rl k+1 + ^ sin{2kft*)/rf} . 

Thus to leading order in 1/n, 

(n - l)/32n+2«2n-2(>*) - (n + l)/32n-2«2n+2 (c* ) ~ -8(n 2 - sin(4^*)/r 4n . 
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Since 7r/4 we can expect that it is a good approximation to choose A such 
that the left hand side vanishes. Doing so, we find a spectacular improvement in 
convergence, cf. fig. 5, which nicely provides further confirmation for our theory. By 
looking at running averages of four points, and assuming no systematic shift from 
neglecting the right hand side of the above, we extract v = .689457(8). Similarly 
we find to = .5955(5). 
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